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1. Introduction 

Yang-Mills theory in the Coulomb gauge has recently drawn a renewed attention, both in the 
continuum [jl], ^ and on the lattice [Q, ||, ^]. In the continuum at least, this interest is mostly due 
to the remarkable fact that GauB' law can be resolved explicitly in Coulomb gauge, which gives 
the remaining vector potential A a very intuitive notion similar to electrodynamics [Q]. Recent 
variational approaches in the Schrodinger picture even support the idea of a constituent gluon [jl], ^j, 
which is almost non-interacting in the infrared and thus completely determined by its dispersion 
relation ffl(p), i.e. the (inverse) equal-time gluon propagator D(p) = jto(p) -1 . 

The obvious drawback of the Coulomb gauge is that manifest Lorentz invariance is lost at inter- 
mediate stages; it may only be recovered at the end of the calculation. Perturbatively, this problem 
is reflected in the (tree-level) propagators of some fundamental fields, which are instantaneous in 
time so that many loop integrands are independent of the temporal loop momentum component 
Icq. Such integrals are notoriously difficult to regulate with conventional techniques, though they 
are believed to cancel in the full theory [^]. Still, the issue of renormalisation in Coulomb gauge 
remains cumbersome, even at the one-loop level 

Similar problems arise on the lattice as well. While initial studies of the gluon propagator 
in Coulomb gauge displayed almost perfect scaling [Q[5|], recent studies using improved gf. tech- 
niques indicate that the quality of gauge fixing has a significant impact on Green functions; in 
particular, substantial scaling violations may result [Q]. The same conclusion has been drawn ear- 
lier in Landau gauge, where careful gauge fixing may alter the infrared behaviour of the propagator 
quantitatively by as much as 20 % [f7|]. 

Even more severe descrepancies arise in the comparision of early lattice results with the vari- 
ational approach mentioned above. While both methods show good agreement in D = 2 + 1, their 
results in D = 3 + 1 differ qualitatively, both in the infra-red and the ulta- violet: 
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All these findings emphasise the need for a thorough corroboration of lattice results in Coulomb 
gauge, in particular with regard to the quality of gauge fixing. In the present talk, I will present 
the first results in this program, viz. the equal time gluon propagator in D = 2 + 1 and D = 3 + 1. 
Further studies on the ghost propagator and the Coulomb form factor are currently underway and 
will be presented elsewhere. 

The plan of this talk is as follows: In the next section, I will briefly discuss our gf. techniques 
and demonstrate that they are effective in reducing the Gribov problem which is at the heart of 
most gf. issues. Section three presents our findings for the gluon propagator. Some of this data 
is still preliminary, and so is the quantitative analysis, but our results so far imply both scaling 
violations in the UV and a significant suppression in the IR. The last point improves the qualitative 
agreement with variational studies, although the quantiative agreement is still unsatisfactory. In the 
last section, I will conclude with a brief summary and outlook. 
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2. Gauge fixing techniques 

Coulomb gauge on the lattice can be denned as the maximisation of the functional 1 

Here, Un (x) are the link variables, the sum over x runs over all sites in a fixed time-slice t = const 
and the maximisation is along the gauge orbit, i.e. with respect to all gauge rotations Cl(x,t) of the 
link field Un(x). As indicated, the Coulomb condition F t = max can be implemented at each time- 
slice t independently. This leaves a residual invariance of space-independent but time dependent 
gauge transformations £l(t), i.e. a global gauge rotation in every time slice. 
For the equal-time gluon propagator 2 

3 3 

D(p) ~ / d 3 xe^<^ £ £ (A? (x,f)Af (j,t)) = Ipr 1 + 0(h) (2.2) 

J i=\ c=\ 

the residual gauge fixing is irrelevant and it is sufficient to fix only the time slice in which the 
measurement is taken. This is no longer true for other correlators such as the Aq — Aq propagator 
related to the static Coulomb potential. Moreover, recent perturbative studies [^] indicate that 
possible scaling violations in D(p) may be attributed to the loss of covariance at equal times; it will 
then be necessary to consider the full gluon propagator at all (unequal) times, and Coulomb gauge 
fixing at all time slices must be augmented by a suitable choice for the residual symmetry. 

The Gribov problem, which is at the heart of most g.f. issues, can be expressed as the fact 



that (gjj) has many local maxima which may, however, give inequivalent contributions to non- 
gauge invariant quantities such as the Green functions. Uniqueness can be enforced by searching 
for the global maximum of ( |2.1[ ), an NP-hard problem. Our strategy to reduce the influence of 
Gribov copies is to prepend the standard (over)relaxation algorithm by an initial preconditioning 
step combined with multiple Gribov repetitions from random starts. This method is a less expensive 
substitute for full simulated annealing and works well for small to medium size volumina up to 
V « 36 4 . 



2.1 Preconditioning 

The periodic boundary conditions on the lattice allow for a somewhat larger symmetry than 
just the periodic local gauge rotations. This is well-known from the SU (2) lattice center symmtry: 
In this case, one multiplies all links Uo(t,x) pointing out of a fixed time-slice 3 t = const by (— 1). 
This construction flips the sign of all Polyakov lines, but it leaves all plaquettes (and thus the action) 
invariant; it is therefore a genuine symmetry of the system. In Landau gauge, one can generalise 
this construction to all four directions, giving a total of 2 4 possible combinations of flips 

'For simplicity, we work exclusively with the colour group G = SU (2). 

2 Gauge potentials are extracted from the link variables in the usual fashion via an 0(a 2 ) improvement of the basic 
formula A M = ^ [ty (x) - (x)] . 

3 The actual location of the time slice / is irrelevant, since a center flip at a different time-slice t' can be decomposed 
into a flip at / followed by a strictly local, periodic gauge transformation. 
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Figure 1: The equal-time gluon propagator at the smallest non-zero lattice momentum, measured as a 
function of the number iV of Gribov repetitions. Data was collected on a 24 4 lattice with j3 = 2.15 (left) and 
j5 = 2.20 (right); a total of 200 thermalised configuration were analysed for each data point. 

In the Coulomb case, the gauge fixing is carried out in a fixed 3D time-slice, i.e. the flips are 
only carried out in spatial directions, and only the 2D sub-planes perpendicular to a given direction 
(at fixed t) are flipped. The preconditioning consist of trying all 2 3 twists to maximise F t [U } prior to 
the acutal relaxation step. This can be viewed as a non-local update representing a large symmetry 
transformation that no local relaxation algorithm is likely to find. Flips can also be interspersed at 
any time during relaxation, although they are most efficient early on, when the algorithm has not 
yet converged onto a target maximum. 4 

2.2 Multiple Gribov repetitions 

The gf. sequence consisting of preconditioning, relaxation and overrelaxation can be repeated 
multiple times with random starting points. This inspects different regions of the search space and 
converges to distinct Gribov copies. What makes this repetition effective is that a relatively small 
number N of copies gives a large increase in the gf. functional, while subsequent repetitions beyond 
a certain plateau point do not give any substantial improvement within reasonable computation 
time. 

This can be seen in figure 1, which plots the equal-time Gluon propagator D(p) at the smallest 
non-zero lattice momentum, as a function of the number ,/V of Gribov repetitions. The net effect of 
the improved gauge fixing is generally to suppress D(p m j n ). Even for ,/V as small as N = 2, ... ,5, 
the corrections are in the range of 10%. Further copies give smaller corrections; it is then a matter 
of experiment to find the optimal tradeoff between CPU time and gf. quality. The optimal ,/V will 
depend quite sensitively on the lattice size and other simulation parameters. In fig. 1 one can see 
the plateau setting in rather quickly, while our largest lattices (V = 36 4 ) required up to ,/V = 30 
repetitions. 

4 The (over)relaxation algorithm is iterated until the local gf. violation, i.e. the (maximal norm at all sites x of the) 
local gradient of (EH!) is smaller than 10~ 13 . 



# Gribov copies 
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3. Results 

3.1 D = 2 + l 

In this case, our findings in rig. 2 are in fair agreement to previous lattice calculations [Q, |5||. 
Our improved gauge fixing scheme has again the tendency to suppress the gluon propagator in the 
infra-red, but since Z>(|p|) — » at small |p| even without gf. improvment, the qualitative behaviour 
of the gluon propagator is unchanged. 

In the UV, we observe scaling in the sense that the various propagator curves for different 
values of the coupling j3 can be multiplied by a momentum-independent factor Z(j8) such that all 
curves coalesce to a single line. There is a tendency for the scaling to be less perfect than without 
the gf. improvement, but this is well below the error bars of our numerical simulation. 

Quantitatively, the suppression of the gluon propagator in the infra-red is as large as 10% - 
15%. To fit the curve in the deep IR and UV region, we have placed two cuts on the data. In the 
IR, a power ansatz yields 

D(|p|) = |p| a -(c 1 +c 2 |p| 2 + ---), a « 0.85(10). (3.1) 

Since the curve flattens towards the maximum, the exponent a is somewhat depending on the exact 
location of the IR cut. At A = 0.5 GeV, we have a = 0.81, while it increases to the above value 
a = 0.85 for A = 0.4 GeV. With our present lattice sizes, we cannot go much lower with the IR 
cut, but the present trend does certainly not rule out the value a = 1 preferred by Hamiltonian 
approaches [|J]. 

In the ultra-violet, a power-law decay 

£>(|p|)Hpr r , 7-1.5(1) (3.2) 

is possible, but the exact value of the exponent y depends quite sensitively on the location of the 
UV cut. A double-logarithmic plot in the deep UV is not a straight line at large momenta, which 
points to sizeable logarithmic corrections. In fact, an ad-hoc ansatz 

D(p) 1 



Pi " m |p| 5 

with 8 fx 0.51 can fit the data equally well. The conclusion is that our present data does not contain 
large enough momenta to distinguish between a logarithmic or a power-like behaviour in the UV. 

3.2 5 = 3+1 

The left panel of fig. 3 shows the results for the largest lattice that we considered. The improved 
gf. scheme is now seen to make a qualitative difference, both in the IR and the UV. 

At low momenta, the propagator is clearly suppressed as compared to less intricate gf. proce- 
dures. The power-law fit explained in the last subsection reveals a IR exponent of 

a » 0.24(12), 

again with significant variations as the IR cut on the data is changed. However, a value a = 0, i.e. a 
gluon propagator going to a constant as p — » [|], |5|] seems much more unlikely than the vanishing 
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Figure 2: The renormalised equal-time gluon propagator for various couplings and lattice sizes. For the 
significance of the two data cuts, see the main text. 

D(0) = predicted by variational calculations [JxT] . On the other hand, the comparision with the 
D = 2+1 case indicates that much smaller momenta must be sampled to rule out one or the other 
option. 

In the UV, the most striking difference to previous lattice results is the absence of perfect 
scaling, i.e. the gluon propagator does not seem to be multiplicatively renormalisable. This can be 
clearly seen in the logarithmic plot in the right panel of fig. 3. In a multiplicatively renormalisable 
situation, we would expect the curves for all couplings /3 to have the same slope at large momenta 
- which is clearly not the case. 

One can now proceed and renormalise anyway such that a common curve can be observed in 
one /^-region or the other (the right panel of fig. 3 has been renormalised to fit well in the IR). In 
particular, one could try to fit the deep UV region, at the expense of sacrificing a common curve in 
the IR. From such a fit, it is even possible to extract a power-like behaviour 

D(fi)~\p\~ a , a = 1.57. 

which is in fair agreement with ref. [Qj. Our present data, however, does not warrant such a pro- 
cedure. In particular, an ad-hoc logarithmic ansatz as in the last subsection would work equally 
well. To summarize, the scaling violations displayed by our improved gf. scheme are so severe 
that any attempt to extract a consistent UV behaviour from a multiplicative renormalisation seems 
ill-adviced. 

Comparable problems with renormalisation were also found in other studies employing im- 
proved gf. schemes. Continuum perturbation theory [^] attributes the scaling violations to the 
instantaneous nature of the propagator considered here. 

4. Summary and conclusions 

In this talk, I have presented first results for the equal-time gluon propagator measured in an 
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Figure 3: Left panel: The equal-time gluon propagator for various values of the coupling constant. The 
gauge fixing includes preconditioning and a minimum of 30 Gribov repetitions for each measurement; mul- 
tiplicative renormalisation focused on the IR data. Right panel: The same data in a logarithmic plot. 



improved Coulomb gauge fixing scheme. The general observation is a significant suppression of 
the propagator in the infrared, and a loss of scaling at very large momenta. Although the numerics 
is not fully compelling, the IR data points to D(0) = as a likely scenario even for D = 3 + 1. The 
failure of multiplicative renormalisation in the UV has also been observed in other studies treating 
Coulomb gauge with improved gf. techniques; in perturbation theory, this failure can presumable 
be attributed to a loss of covariance for the equal-time propagator. 

To make the present numbers more convincing, we have to go to smaller momenta, which may 
involve a simulated annealing step in the gf. pipeline. To get a handle on the scaling issue, it would 
also be interesting to study the gluon propagator at non-equal times, using a complete gauge fixing 
that also destroys the residual symmetry in Coulomb gauge. Further investigations involve the 
ghost propagator and the Coulomb form factor, which are of immediate relevance for the physics 
of the gauge system. These studies are currently underway and will be presented elsewhere. 
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